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Abstract 

The photon radiation in the initial state lowers the energy available for the e + e~ 
collisions; this effect is particularly important at LEP2 energies (above the mass 
of the Z boson). Being aligned to the beam direction, such initial state radiation 
is mostly undetected. This article describes the procedure used by the DELPHI 
experiment at LEP to estimate the effective centre-of-mass energy in hadronic events 
collected at energies above the Z peak. Typical resolutions ranging from 2 to 3 GeV 
on the effective center-of-mass energy are achieved, depending on the event topology. 
PACS codes: 07.05. K, 13.10.+q, 13.65.+i. Keywords: Data Analysis, Photon Ra- 
diation, Initial State Radiation. 
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1 Introduction 



The increase in the LEP beam energy above the Z in 1995 opened a new window 
of higher centre-of-mass (cm.) energies to study e + e~ annihilations. The Initial State 
Radiation reduces the energy available for the physical processes. 

The photons from Initial State Radiation (hereafter called ISR photons) tend to move 
along the beam (z) direction; thus, they are mostly undetected. Figure [I] shows the energy 
versus polar angle distribution of the ISR photon(s) in simulated hadronic events at a cm. 
energy of 161 GeV. The initial qq state was generated using the PYTHIA 5.4 |l| Parton 
Shower Monte Carlo. The polar angle, 8, peaks at values close to zero, and the photon 
energies cluster at E lR — (s — M|)/(2y / s), where s = (2E b ) 2 , E b is the beam energy, and 
Mz is the Z mass. For events with such a hard ISR, the energy in the cm. of the e + e~ 
collision is reduced to the Z mass ("Z radiative return" events). 

This note explains the methods used by the DELPHI experiment at LEP to compute 
the effective centre of mass energy, \fs' \ in e + e~ annihilations into hadrons, by measuring 
or estimating the energy and momentum of ISR photon(s). Two algorithms are used: 

• the first one, based on the jet directions only, is more robust; 

• the second, based on a constrained fit, is more precise. 

The procedures described in this note were used in many papers from DELPHI (see 
for example || §). The DELPHI detector is described in |§, and its performance in p. 

2 The methods 

One looks first for detected ISR photons. Showers unassociated to charged particles 
are searched for in the calorimeters, and they are considered as candidate ISR photons, 
if they fulfill the following criteria: 

• energy deposition larger than 10 GeV; 

• isolation angle larger than 0.3 radians with respect to any charged particle with 
momentum larger than 1 GeV; 

• when the shower energy is smaller than 70 % of the photon energy corresponding to 
a radiative return to the Z, isolation angle with respect to any charged particle (with 
momentum larger than 1 GeV) larger than 0.6 radians and smaller than 1.5 radians 
(1.8 radians if the shower is contained in the small angle calorimeter, the STIC). 

Showers are grouped together, if they are closer than 5 degrees in space. The associa- 
tion is done by adding the energies and assuming as direction the direction of the vector 
sum of the momenta. 

2.1 Method 1: determination of the energy of ISR photon(s) 
from the direction of jets 

The first method (called SPRIME inside the collaboration) uses only the direction of 
jets and of the detected photon(s) - if any. 
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All the particles in the event, except the one(s) considered as ISR photon(s) using 
the criteria described above, are forced to be clustered in two jets, using the DURHAM 
algorithm [[J. The particles are assumed to be massless. 

The method applies energy and momentum conservation to estimate the energy of 
possible ISR photons lost along a direction close to the beam, after having considered the 
candidate one(s) directly detected in the electromagnetic calorimeters. 

If no ISR photons are detected, an event plane containing the beam (z) axis and 
maximizing the sum of the moduli of the projections of the jet momenta is determined. 
A hypothetical photon is assumed to be along z, in the positive or negative direction, 
respectively if the sum of the jet polar angles is larger or smaller than 180 degrees. The 
directions of the two jets, projected onto the event plane, are used to constrain the photon 
energy. 

If only one ISR photon is measured in the detector and the event is planar, i.e., if 
the sum of the internal angles between the jets and the photon exceeds 345 degrees, 
the photon energy is calculated by energy-momentum conservation in the event plane 
containing the photon direction and maximizing the sum of the moduli of the projections 
of the jet momenta. For aplanar events at least one particle has been lost. The algorithm 
assumes that at another photon has been radiated inside the beam pipe and makes use 
of the energies of the jets and of the measured ISR photon to constrain the hypothetical 
photon energy. 

If more than one ISR photon is measured, energy-momentum conservation is applied in 
the event plane to determine the energy of an unseen photon along z and the momenta of 
the two jets, using as an input the direction of the jets and the energies and the directions 
of the detected photons. 

2.1.1 Results 

The estimated effective cm. energy (V?), for hadronic events (selected using the standard 
DELPHI hadronic criteria, see (?]) collected from the run around 188.7 GeV during 1998, 
is shown in the upper plot of figure ^| (data points) and is compared to a sample simulated 
with PYTHIA and followed through the complete simulation of the detector, DELSIM ||. 
The results show a satisfactory description of the data by the simulation. 

High efficiencies and purities are obtained for selecting "high energy" events (defined 
as those with effective energy above 0.9-y/s), and for selecting Z radiative return events. In 
the following, we shall call "low energy" events those with effective energy below 0.9^. 

The generated effective cm. energy for simulated events at beam energy of 80.5 GeV, 
selected requiring a SPRIME result above 0.9^/i, is shown in figure [3] (hollow histogram). 
The hatched histogram shows the same distribution for a subsample of those events in 
which the energy seen in the detector is above M z + (y/s — Mz)j1 = 126 GeV. With this 
cut the small contamination due to events with effective cm. energy near the Z mass is 
drastically reduced. 

The efficiency to select high energy events was computed to be 88%, with a purity of 
the selected sample of 92%. When a larger purity is required, a cut on the total energy 
of the event may be applied. For the cut illustrated in the previous paragraph the purity 
was increased to 96%, with an efficiency of 74%. 

Using simulated events with PYTHIA at cm. energies of 130 GeV, 161 GeV and 
188 GeV respectively, the efficiencies and purities for selecting "high energy" events, as 
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a function of the cut value in \fs' / y/s, are shown in figure f|. 

Figure || shows the distribution of the low energy events, simulated at cm. energy of 
188 GeV, as a function of the difference between the SPRIME estimate and the generated 
effective energy. The Full Width at Half the Maximum (FWHM) of this distribution 
divided by 2.35 is used as an estimator of the resolution of SPRIME, and is shown on the 
left plot of figure |B| as a function of the DURHAM y cut value needed to distinguish 2-jet 
from 3-jet events. The right plot of figure || shows the average offset of the method, defined 
as the peak value of the distribution of the difference between the SPRIME estimate and 
the generated effective energy. 

2.2 Method 2: determination of the energy of ISR photon(s) 
from a constrained fit 

The second method (referred to inside the collaboration as SPRIME+) uses a con- 
strained fit taking into account the 4-momenta of the measured jets and photons and 
assuming one additional photon inside the beampipe. This method is able to include any 
number of jets and photons and is therefore more appropriate than the previous method 
in events with more than two jets. The fit also returns an estimated error on the effective 
cm. energy, and a x 2 indicating the goodness of the fit. 

All particles that have not been identified as ISR photons are clustered into a "natural" 
number of jets using the DURHAM jet algorithm with a y cut fixed to 0.002. This y cnt 
value was chosen to optimize the energy resolution; at this value, PYTHIA reproduces 
the data well. 

Then two kinematical constrained fits are performed, corresponding to the hypotheses 
that an undetected photon went along the positive or negative z-direction (along the 
beam). There are two constraints from the transverse momentum balance, J^Px = and 
Y^Py — 0, and one from conservation of the total energy and of the momentum along the 
beam direction: 

Y,E = y/s-Y,Pz or Y,E = y/s + J2Pz 

assuming J^Pz > or J^Pz < respectively. 

The fitted jet momentum, py , was projected onto a set of axes with one component 
parallel to the measured jet momentum, p~j m , and two transverse components. The parallel 
component was described by a rescaling factor, exp(aj), while the transverse components 
were described by parameters multiplying perpendicular momenta fixed to 1 GeV/c: 

Pj f = e Maj)fj m + bjPj b + Cjpf. 

To determine p~j b and p~j C , the jet energy, Ej m , and thereby also to a good approximation 
the jet mass, was rescaled with the same factor exp(aj) as the jet momentum; the fitting 
algorithm then minimized a \ 2 '- 

while forcing the fitted event to obey the constraints. 
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The errors were parametrized as a function of the polar angle 6j of the jet with respect 
to the beam: 

a = 0.15 + 0.40 • cos 4 ^- 
a a . = 0.15 + 0.40 • cos 4 ^- 
cr b . = a Cj = 1.8 + 1.08 ■ cos 4 ^- GeV 

Photons are treated in the same way as jets, but with smaller estimated measurement 
errors: 

a = 
a aj = 0.04 
a bj = a Cj = 0.01 GeV. 

When, after the fit, J^Pz has the same sign as assumed in the constraints, the solution 
is said to be "physical" , otherwise it is called "unphysical" . The solution is chosen that has 
the smallest x 2 an d is physical. If both solutions are unphysical the measured energy and 
momentum balance indicate that there was no significant photon lost inside the beampipe. 
In that case a new constrained fit is performed, using the four constraints of momentum 
and energy conservation with J^Pz fixed to zero. 

The fitting algorithm used is similar to the method described in detail in |§. 

In case the standard fit as described above does not converge, a special procedure 
is applied to try to recuperate the event. The longitudinal errors are relaxed for the 
measured jets: 

a = 0.75 + 2.0 • cos 4 ^- 
a a . = 0.75 + 2.0 • cos 4 ^- 

and for the photons: 

a = 
a a] = 0.057. 

Then a series of fits is done with the 4 constraints of energy and momentum conservation 
after assuming an additional photon along z, with momentum p z = -100, -95, -90 ... +95, 
+ 100 GeV along the beam. 

The fit with the smallest \ 2 is chosen. If there are at least three fits around the mini- 
mum that converged, interpolation with a parabola is used to improve the determination 
of the p z with the lowest \ 2 - A final fit is done with p z constrained to this value. 

2.2.1 Results 

The performance of the constrained fit method was studied using PYTHIA qq7 simulation 
at cm. energies of 130 GeV, 161 GeV and 188 GeV, including the full detector simulation 
DELSIM. 

The constrained fit turns out to be quite robust. The standard fit does not converge 
in 0.4% of the events, in which case the recuperation method is used. This method is able 
to recuperate 60% of the events that failed. 
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The resolution obtained with SPRIME+ can be seen from Figure [5]. The improvement 
in resolution and offset with respect to SPRIME is evident especially for events with a 
high y cut from 3 to 2 jets, for which the assumption of 2 massless jets is no longer valid 
(at a cm. energy of 161 GeV, 46% of the events have y cut > 10~ 3 ). This is shown for the 
reconstruction of Z radiative return events in figure || 

The improvement in resolution leads to higher efficiencies for selecting events with 
y/s* larger than 0.9a/s (figure |4j). For very high purities, however, there is hardly any 
improvement. This is due to a side effect of the improved resolution for 3-jet like configu- 
rations: when an energetic ISR photon is radiated inside the detector acceptance but not 
tagged as an ISR photon, it is be treated as a third jet in SPRIME+, giving an estimated 
effective cm. energy larger than 0.9^/i. As SPRIME has a worse resolution for 3-jet like 
events, a significantly larger part of events of this type of events miss the cut, resulting 
in a higher purity. 

With a cut requiring SPRIME+ larger than 0.9-y/s, the efficiency and purity for events 
with generated -v/s 7 larger 0.9-^/i for the PYTHIA MC sample at 161 GeV are 94% and 
90% respectively. 

In figure |2| the data taken at 189 GeV is compared to simulation. The high energy 
peak is well seperated from the radiative return peak. The figure also shows the improved 
resolution for the background which mainly consists of W pair events. The multi-jet 
events coming from the fully hadronic decay are reconstructed very well at the full cm. 
energy. The large tail at low ^Js consists of hadronic-leptonic and fully leptonic W pair 
events containing neutrinos that are not included in the reconstruction hypothesis. 

3 Summary 

The methods used by DELPHI for the determination of the effective center-of-mass energy 
in hadronic events at LEP 2 were presented. 

The "high energy" peak is very well separated from the radiative return to the Z, 
allowing for efficiencies and purities well above 80% to be obtained for selecting events 
with effective energy greater than 0.9a/s ("high energy" events) or in the region of the Z 
mass ("Z radiative return" events). 

Typical resolutions ranging from 2 to 3 GeV on the effective center-of-mass energy are 
achieved, depending on the event topology. 
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Figure 1: Energy as a function of polar angle for ISR photons at cm. 
energy of 161 GeV. Photons below 2 degrees are not plotted. 
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Figure 2: SPRIME (top) and SPRIME+ (bottom) results for the data collected during 
early 1998 at a cm. energy around 188 GeV (black bullets), compared to the simu- 
lated events at the same energy using PYTHIA. The data correspond to an integrated 
luminosity of about 54 pb _1 . 
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Figure 3: Generated effective centre of mass energy for simulated 
events with a SPRIME result above 0.9-^/s (solid line), and the same 
distribution for events with a SPRIME result above 0.9^ and requiring 
that the total energy seen in the event is above 126 GeV. 
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Figure 4: Efficiency (black bullets) and purity (white squares) for 
selecting events with generated effective centre of mass energy above 
0.9a/s as a function of the cut applied on the ratio of the SPRIME 
(left) or SPRIME+ (right) result to the LEP cm. energy. The results 
are shown for simulated events at cm. energies of 130 GeV (two upmost 
plots), 161 GeV (two central plots) and 188 GeV (two lowmost plots). 



10 




Figure 5: Solid line: distribution of PYTHIA simulated events at an 
energy in the centre of mass of 188 GeV as a function of the difference 
between the SPRIME result and the generated effective cm. energy. 
Dashed line: same for SPRIME+. 




Figure 6: Left: Resolution of SPRIME and SPRIME+, as a function of 
the minimal y cut needed for DURHAM to reconstruct two jets. Right: 
Average offset for SPRIME and SPRIME+, as a function of y cut . 
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